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Levitated nanospheres in optical cavities open a novel route to study many-body systems out of 
solution and highly isolated from the environment. We show that properly tuned optical parameters 
allow for the study of the non-equilibrium dynamics of composite nano-particles with non-isotropic 
optical friction. We find optically induced ordering and nematic transitions with non-equilibrium 
analogs to liquid crystal phases for ensembles of dimers. 



Introduction - The interaction between light and mat- 
ter has been one of the central driving forces behind 
recent developments in condensed matter physics with 
nanoparticles [lj. Optical tweezers [2j|3] and confocal mi- 
croscopy [4] have made it possible to study many-body 
systems of nanoparticles in solutions in real-time and 
with single-particle resolution [5-8 . Recently, it has been 
proposed to optically levitate and cool single nanospheres 
inside an optical cavity [9UT2] . While in the realm of soft- 
matter physics, this may provide an alternative to con- 
finement in a solution [T3] . from a quantum-optics point 
of view such a set-up provides a versatile alternative to 
conventional optomechanical systems [T4UT6] . Combined 
with optomechanical cooling and trapping techniques of 
single particles, this may even open the possibility to 
study fundamental aspects of quantum mechanics with 
mesoscopic objects [TTJ [T71 [18]. Here, we focus on the 
dynamics of many, interacting particles in the presence of 
optomechanical cooling. While many-body systems with 
non-uniform cooling have been studied with atoms and 
ions [19] , the possibility to create complex structures with 
nanospheres offers completely new opportunities to study 
pattern formation and self assembly. With novel synthe- 
sis methods it is now possible to design compound struc- 
tures ranging from dimers to networks of nanospheres 
connected by spring-like biomolecules [20-26] . A distinc- 
tive feature of the self-assembly of composite particles 
is that the emerging patterns are characterised not only 
by their positions, but also by their individual orienta- 
tions [27l|28]. The non-equilibrium self-assembly of such 
nano-structures in the presence of non-isotropic optical 
cooling is an open question, and holds the promise of new 
means to optical control of pattern formation and novel 
non-equilibrium liquid crystal phases. 

In this Letter, we study the dynamics and self- 
assembly of levitated nanosphere-dimers in presence of 
optomechanical friction inside a two-mirror cavity. The 
particles are subject to thermal forces and coupled to 
a cavity mode which is driven by an external laser and 
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damped by the cavity-decay. The optomechanical inter- 
action gives rise to an optical potential and a cooling force 
along the cavity axis. By compensating the potential 
with a second optical mode, the remaining optomechani- 
cal effect is friction along the cavity axis. Fig.[l]illustrates 
the system we have in mind. The nanospheres are har- 
monically trapped and confined to the xy plane, where x 
is the cavity axis. We show that the steady-state orienta- 
tion of a single dimer is non-uniform in presence of optical 
friction. Remarkably, the full range of preferred orienta- 
tions from to 7r/2 is accessible by appropriately tuning 
the experimental parameters. In a many-body system, 
the presence of additional direct interactions between the 
individual nanospheres, leads to competition between the 
natural triangular ordering of a two-dimensional crystal 
and dissipation-induced ordering. Compared to other ap- 




FIG. 1: (a) Levitated dimers composed of dielectric nano- 
spheres in a laser driven optical two-mode cavity with laser 
strengths Qi,2, a decay rate k and an additional harmonic trap 
Vtrap. The orientation of the dimers is denoted 0, their equi- 
librium separation xq and their frequency v. The direct pair 
interaction between the nanospheres leads to liquid crystal 
phases with rich dimer patterns, (b) The particles are in con- 
tact with a thermal bath with coupling 7. The cavity-particle 
interaction results in strong additional dissipation in one spa- 
tial direction only, which acts as an effective zero tempera- 
ture bath 7 pt- ( c ) Comparison of the relevant time scales: 
We assume that 7 op t is faster than 7 and find 3 qualitatively 
different regimes, which give rise to different orientations, de- 
pending on the relative timescales. 
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proaches to dynamical ordering of dimers, such as shear 
[29] or static electric fields [30], this method offers ad- 
ditional advantages: (a) the orientation depends on the 
frequency of the vibrational mode of the dimers which al- 
lows for individual ordering in multi-species systems and 
(b) in addition to ordering at the level of single parti- 
cles non-uniform friction also leads to novel liquid crystal 
phases at the many-body level (see Fig. [3|. We identify 
three relevant time scales in this system: the frequency 
of the vibrational mode of the dimers the optomechan- 
ical damping rate 7 opt and the rate of thermalization 7. 
For limiting cases we present analytic results on the or- 
dering of individual dimers and numerical results on the 
non-equilibrium many-body dynamics. 

Model - We consider a system of N/ 2 dimers consisting 
of N nanospheres, trapped inside and optical cavity. The 
Hamiltonian is decomposed as H = H sys + H om with 



(i) 



^^-(l X 2j-l - X 2j| -%0) 2 +r ^^nt(|Xi -Xjl) 



J = l 



the system Hamiltonian. Here, m is the mass of the 
nanospheres, Xj and pj are the position and momentum 
of particle j, respectively, Vtrap is the trapping potential, 
v the frequency of the vibrational mode of the dimers, Xq 
is the equilibrium separation and VJ nt is the direct dipolar 
pair interaction, which can be tuned by the parameter 
To [6]. In a frame, rotating with the laser drive, the 
optomechanical Hamiltonian is given by 



H om = -A(xi ... xLj)\a\ 2 + Sl(a + a*) + V c 



(2) 



where a(t) [v Js] is a normal variable, which describes 
the dynamics of the optical mode, Q [y/ J/s] character- 
izes the drive strength, V c is a second compensating po- 
tential, and A = uod — gj(xi ... Xj) is the detuning of 
the laser drive from the optical resonance. The inter- 
action derives from the electric polarizability a p of the 
nanospheres [9j[l0], which, for sub- wavelength particles, 
gives rise to a position-depended cavity resonance fre- 
quency cj(xi ... Xj) = uj - (gC/2) l F (xj)| 2 , where uj 
is the bare cavity frequency. F(x) is the normalized mode 
function of the cavity mode, C is the mode volume and 
g = a p uo/C is the optomechanical coupling strength, 
which is proportional to the volume of the particle. We 
assume that the cavity mode can be approximated by 
a standing wave F(x) c± e^/2/C sin(fcx), where k is the 
wave number and e is the polarization in the yz plane, 
so that A(xi ... Xj) = A(xi ... Xj). 

The full dynamics of the optical part, including decay 
from the cavity, is described by the equations of motion 
for the optical amplitude a = (iA(xi ... Xj) — n)a + iVt, 
where k is the field decay rate. The drive can be elim- 
inated by the transformation a — » a + a' with a = 
fi/(A + in) the average amplitude and A = (uod — uo) 



the detuning from the bare cavity resonance. The equa- 
tion of motion for the transformed optical amplitude is 
a! — (iA — k)oJ + i{A{x\ ... xj) — A) (a + a'). Sub- 
stitution of a — » a + a' in the first term in Eq. ([2| 
gives rise to four terms. One of them corresponds 
to an effective optical potential along the cavity axis 
V{x\ ... Xj) = —\a\ 2 A(x\ ... Xj), which can be compen- 
sated by V c of the second optical mode. For the choice 
of a mode separated by one free spectral range from cj, 
the mode function in the focal range of the cavity is ap- 
proximately given by G(x) ~ r)y/2/C cos(kx), where 77 is 
the polarization. When e and rj are orthogonal, the two 
modes do not interfere and the effective potentials add 
up, i.e. V(xj) — — g(\a s \ 2 sin 2 (kxj) + |a c | 2 cos 2 (kxj)). 
For \a s \ = |a c |, the potential is independent of Xj and 
the forces cancel. The second mode is driven on reso- 
nance and therefore does not lead to additional damping 
(or amplification). 

We find the equations of motion of the nanosphere po- 
sitions and momenta 
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dA(xi ... xj) 

dxi 



. (3) 



The cavity decay rate n sets a finite time scale for the 
cavity to respond to changes in the particle positions so 
that the optomechanical feedback does not only depend 
on the particle positions, but also on the particle mo- 
menta along the cavity axis [14]. This gives rise to am- 
plification (for A > 0) or damping (for A < 0), respec- 
tively. Optomechanical cooling of nanospheres has been 
studied previously in the Lamb-Dicke regime [9 , while 
here, we focus on a different regime, similar to Ref. [33] , 
in which the particles move almost freely and the rele- 
vant frequency stems from the modulation of the optical 
feedback force at co> mo d = 2kp/m. The damping force has 
resonances for A = ±|cj mo d|, while in the regime we are 
interested in, \p\ <C mA/(2fc), the cooling rate is approx- 
imately constant 



7opt 
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(4) 



Finally, the equations of motion take the form of a mod- 
ified Langevin equations with 



dH< 
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sys 



dy 



■ ra(7 + 7o P t)x + 6c 



(5) 



Here, 7 is the rate of thermalization and £ x and £ y are 
mutually uncorrelated Langevin forces, as characterized 
by = k B TmjS(t - t'). 

Results - We first focus on the dynamics of a single 
dimer described by Eq. Separating the dynam- 

ics into the trivial center-of-mass motion and the rel- 
ative coordinates (xi — X2)/V2 = (x,y), and to the 
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FIG. 2: Steady-state probabilities of the orientation as a 
function of the dimer frequency from numerical integration of 
the second order Langevin equation Eqs. ([5| for cooling rate 
ratios 7o P t/7 = 10 (a), 7o P t/7 = 100 (b), and 7o P t/7 = 1000 
(c). The comparison of the analytic expressions Eq. (T6j) for 
small frequencies (d dashed lines) and numerical results (d full 
lines) show good agreement over a wide range of frequencies, 
(e) Analytic results from Eq. ^ for a range of frequencies 
with 7o P t/7 = 100 predicts the transition from = orien- 
tation to — 7r/2. For larger frequencies the model breaks 
down as it predicts = tv/2 for large v. (e). In the rigid-rotor 
regime, the analytic approximation Eq. ([8} is in agreement 
with the numerical simulation as shown for various cooling 
rates r c = 7o P t/7 (f)- 



extent that optomechanical coupling between the par- 
ticles can be neglected, the non-linear force is dI ^ a = 
— rmv 2 (l — (2xq/(x 2 + y 2 )) 1 / 2 ) , where r = (x,y) and 
xq is the dimer separation. Fig. [2] (a-c) show the steady 
state solutions of the orientation of single dimers as a 
function of the dimer frequency for various cooling rates 
from numerical integration of Eq. (J5|. Remarkably, when 
exposed to uni-directional friction, a loosely connected 
dimer as well as a rigid rotor tends to align orthogonal 
to the direction of friction, whereas a dimer of moder- 
ate stiffness aligns parallel to it. This can be understood 
from three competing effects which derive from the order 
of the relevant time scales in the system (see Fig. [lj): 
7 _1 , 7^ and v~ x . While 7 and 7 opt set the scales of ther- 
malization and non-uniform friction, v sets the time scale 
at which the degrees of freedom mix due to the non-linear 
nature of the force term. Assuming that 7 op t > 7, there 
are three limiting parameter regimes: (i) 7 opt ^> 7 ^> v, 
(ii) 7o P t > v > 7 and (iii) v > 7 > 7 opt . In the follow- 
ing we analytically study these limiting cases and give an 
intuitive explanation of this remarkable non-equilibrium 
ordering phenomenon. 

In (i) and (ii), 7 opt is the largest scale. Due to the 
resulting separation of time scales of the motion in the 
x and y directions, the steady-state distribution is of 
the form P(x,y) = P(y\x)P(x). Here, P(y\x) is the 



normalized steady-state distribution of the fast y direc- 
tion given a fixed x. Consequently, the average energy 
(V)x — ST OQ V( x ->y)P(y\ x ) determines the distribution 
P{x) for the x direction via a Markov process with the 
scaled temperature T x — Tj/(j + 7 op t), thus 



P(x,y) = 



-V(x,y)/(k B T) 



-(V) x /(k B T x ) 



f°° dv e~ v ( x >y)/( k BT) f°° dx e -< y >-/(^ T -) 
J— 00 y J— 00 



(6) 

In regime (i), the spatial fluctuations are much larger 
than #o, so that V(x,y) ~ mv 2 {x 2 + y 2 )/2 and the in- 
tegrals can be evaluated analytically. The distribution 
P(4>) = Jq°° rdrP(r cos 0, r sin 0) for the orientation of 
the dimer 6 is found as 



p{4>) 



\/7(7 + 7o P t) 



27r((7 + 7o pt ) cos 2 <j> + 7 sin <j>) 



(7) 



with the maximum for = 0. In case (ii), the harmonic 
approximation of V(x,y) breaks down, and Eq. ([6l is 
evaluated numerically, with the results shown in Fig]2|d 
and e) which are in agreement with numerical integration 
of Eq. ([5} shown in Figj2jb). 

In the rigid-rotor regime (iii), motion in the radial di- 
rection is suppressed, so that ^x 2 (t) + y 2 (t) = xq and 
we can derive a Langevin equation for <j> alone: mxo(j) = 
-772,(7 + 7 op t cos 2 <t>)<j) + ^ where ^ = £ x sin </> + £ y cos <j> 
so that {£<t>(t)£><f>(t')) = kBTm^S(t — t') is independent 
of (j). This equation describes thermal motion of the 
orientation of the dimer with angle- dependent damping. 
Since there is no conservative force, the motion is over- 
damped, so that <j) ~ 7(0)0, where 7(0) = 7+7 op t cos 2 0, 
and the Fokker-Planck equation reduces to dP/dt = 
(7 + 7opt c °s 2 (j))~ 1 dP/d(j). The stationary solution is 



p(0) = 



7 + 7o P t cos 2 <p 

27T7 + 7T7op t 



(8) 



and has a maximum at = 0. 

This motivates the following intuitive picture: The 
dominant mechanism in (i) is purely geometrical. Here, 
the distribution P(x, y) is a Gaussian that is squeezed 
in the direction of cooling x and, therefore, the most 
likely orientation is = 0. In the rigid-rotor case (iii), 
in which no orientation-dependent energies are involved, 
the anisotropy of the steady-state orientation is due to a 
purely dynamical effect. In this case, the dimer is dynam- 
ically attracted to orientations for which the fluctuations 
in the angular direction are suppressed. The opposite ori- 
entation is reached in the intermediate regime (ii), when 
a third, purely energetical effect, is dominant. The mo- 
tion is mostly confined to configurations of constant x, 
for which the potential energy changes from a double har- 
monic well at x — to a single anharmonic well (V oc x A ) 
at x = ±xq. With the average kinetic energy (T) con- 
stant and using the virial theorem (V) = 2/n(T), with 
n the power of the external potential, we find that the 
probability exp[— /3 X (V}] is largest for = ir/2. 
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Let us consider the experimental feasibility in a con- 
figuration as shown in Fig. [I] The confinement of the 
dumbbells to the xy plane can be provided by an ex- 
ternal standing wave optical trap crossing the Fabry- 
Perot cavity in z-direction. The three relevant time 
scales can be controlled over a large range of parame- 
ters: Dimer frequencies up to v ~ 2tt x 1 kHz can be 
reached e.g. with spring constant k ~ 0.2 pN//im of 
DNA and silica nanospheres with a radius of r = 50 nm 
(mass m = 1.2 x 10 -18 kg). We note that spring con- 
stants smaller by orders of magnitude are possible above 
the persistence length of 50nm [32]. The optical damp- 
ing is provided by a cavity with length L ~ 10 ~ 2 m 
and mode waist of w ~ 10 -4 m. We find the optome- 
chanical coupling g ~ 2ir x 10 4 Hz via the mode vol- 
ume C = (tt/4)Lw 2 . We further assume a cavity finesse 
of T ~ 10 5 , so that k ~ 2 x 10 5 Hz at a wavelength 
of A = 1064 nm. When we further choose A > 5ft, 
a power for the cooling laser of Pdrive = 3 x 10 -5 W 
results in a cooling rate of 7 opt < 2tt x 4 Hz. For 
the thermal environment we assume room temperature 
T = 293 K and 7 = 0.05 Hz, which corresponds for the 
chosen nanospheres to an air pressure of approximately 
10 -5 mbar. Therefore, rates of up to 7 op t/7 ~ 100 are 
possible. Note that even higher rates have recently been 
achieved experimentally by optical feedback cooling [31] . 

Many-body phases - Liquid crystal phases of dimers 
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FIG. 3: (a) Ensemble of dimers with frequency v — 2OO7 (in- 
termediate regime (ii)) in thermal equilibrium are uniformly 
orientated (inset a), (b) The non-isotropic optical friction in- 
duces orientational ordering of individual dimers (inset b) and 
in addition, many-body nematic ordering along the direction 
of cooling, (c) The non-equilibrium phase analog to a liquid 
crystal phase characterized by the pair-correlation function 
g(x). The rigid rotors align orthogonal to cooling with addi- 
tional nematic order in direction of cooling (black) and liquid 
order in ^/-direction (red), (d) Mixture of dimers self assem- 
bles into a crystal with dimer-orientations according to the 
frequencies v. 



have been previously studied in equilibrium and non- 
equilibrium [28] . Non-isotropic friction and the resulting 
ordering may offer novel tools to guide the self-assembly 
towards preferred structures and to study novel nematic 
phases. We numerically study the system described by 
the Eq. with the experimental parameters as given 
above. Fig. [3]^a) depicts an ensemble of dimers without 
optical friction. For this choice of parameters, the sys- 
tem is in the liquid phase and the dimer orientations are 
distributed uniformly. Additional non-isotropic cooling 
with 7o P t/7 = 100 induces a non-equilibrium transition 
to a phase characterized by the single- and many-dimer 
order parameters respectively shown in Fig. [3] (b) and 
(c). The single dimers are still aligned. In addition, the 
interplay between orientation and many-body dynamics 
leads to a remarkable phase with liquid order in the y- 
direction and solid order in the x-direction. With the 
frequency v' = b00v and keeping all other parameters 
fixed, a different pattern with all individual dimers reori- 
entated orthogonal to the direction of cooling (Fig. [3^c) ) 
is found. Again the many-body dynamics leads to order- 
ing along the x-direction, as measured by the directional 
pair-correlation function g x , y (x) = (S((xi — Xj) — x)) de- 
picted in Fig. [3|c). The patterns in Fig. [3] (b) and (c) 
are non-equilibrium analogs of liquid crystal phases in 
equilibrium. For large Tq above the melting temperature 
[6], two-dimensional dipolar particles self-assemble into 
a triangular lattice. The guided orientation of individ- 
ual dimers induced by non-isotropic friction is competing 
with the orientation in the triangular lattice in which the 
dimer orientations are random integer multiples of tt/3. 
This frequency-dependent ordering can be used to cre- 
ate patterns of dimers with different orientations. Fig. 
|3jd) depicts a mixture of dimers with v\ = 500^ and 
with V2 = v. In the absence of direct interactions these 
species order with almost orthogonal relative alignment. 
Increasing the interactions To > r me it above the critical 
melting temperature, the system forms a crystal with tri- 
angular ordering and separate orientational order of the 
two individual dimer species. 

In summary, we have presented a realistic optomechan- 
ical setup that allows for novel dissipative control over the 
orientation of dimers composed of nanospheres. We have 
shown that this approach can be used to prepare non- 
equilibrium analogs of liquid crystals and to study tran- 
sitions in mixtures of multiple species of dimers. The only 
relevant parameters that determine the non-equilibrium 
ordering are the time scales of the vibrational mode of 
the dimers and the rates of thermalization and non- 
isotropic friction. Complex structures of DNA-connected 
nanospheres are of growing interest, and we hope that the 
approach discussed here will provide useful means to op- 
tical control of such systems, complementary to direct 
optical manipulation via optical tweezers. The presented 
mechanism does not rely on any specific properties of 
nanospheres and applies, at least in principle, in general 
to complex structures of dielectric objects, such as viruses 
and bacteria [34]. We speculate that, in the longer run, 
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it may also be applied to complex molecules, and may 
even prove fruitful as a microseeding technique for the 
nucleation of complex molecules. 
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